Fermion-parity anomaly of the critical supercurrent in the quantum spin-Hall effect 
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The helical edge state of a quantum spin-Hall insulator can carry a supercurrent in equilibrium 
between two superconducting electrodes (separation L, coherence length 5). We calculate the max- 
imum (critical) current Ic that can flow without dissipation along a single edge, going beyond the 
short-junction restriction L ^ ^ of earlier work, and find a dependence on the fermion parity of 
the ground state when L becomes larger than ^. Fermion-parity conservation doubles the critical 
current in the low-temperature, long-junction limit, while for a short junction 7c is the same with 
or without parity constraints. This provides a phase-insensitive, DC signature of the 47r-periodic 
Josephson effect. 

PACS numbers: 74.45. -fc, 71.10.Pm, 74.78.Fk, 74.78.Na 



The quantum Hall effect and quantum spin-Hall effect 
both refer to a two-dimensional semiconductor with an 
insulating bulk and a conducting edge, and both exhibit a 
quantized electrical conductance between two metal elec- 
trodes. If the electrodes are superconducting, a current 
can flow in equilibrium, induced by a magnetic flux with- 
out any applied voltage. In the quantum Hall effect, the 
edge states are chiral (propagating in a single direction 
only) and two opposite edges are needed to carry a su- 
percurrent [TH3]. Graphene is an ideal system to study 
this interplay of the Josephson effect and the quantum 
Hall effect in a strong magnetic field [IHS]. 

The interplay of the Josephson effect and the quan- 
tum spin-Hall effect, in zero magnetic field, has not yet 
been demonstrated experimentally but promises to be 
strikingly different [7]. The quantum spin-Hall insulator 
has helical edge states (propagating in both directions) 
that can carry a supercurrent along a single edge. The 
edge state couples a pair of Majorana zero-modes, allow- 
ing for the transmission of unpaired electrons with h/e 
rather than h/2e periodic dependence on the magnetic 

flux 8/9;. 

An h/e flux periodicity corresponds to a 47r-periodicity 
in terms of the superconducting phase difference (j), 
which means that the current-phase relationship has two 
branches I± {(p) and the system switches from one branch 
to the other when (jj is advanced by 2it at fixed total num- 
ber Af of electrons in the system. This is referred to as a 
fermion-parity anomaly, because the two branches have 
different parity cr = ± of the number of electrons in the 
superconducting ground state [8j. 

Josephson junctions come in two types [lOi, depending 
on whether the separation L of the superconducting elec- 
trodes is small or large compared to the coherence length 
^ = hv/A, or equivalently, whether the superconducting 
gap A is small or large compared to the Thouless energy 
Et ^ hv/L. Existing literature [7H51 fTTHT5] has focused 
on the short-junction regime L ^. The supercurrent is 
then determined entirely by the phase dependence of a 
small number of Andreev levels in the gap, just one per 
transverse mode. The phase dependence of the contin- 



uous spectrum above the gap can be neglected. As the 
ratio L/^ increases, the Andreev levels proliferate and 
also the continuous spectrum starts to contribute to the 
supercurrent. Since a is switched by changing the occu- 
pation of a single level, one might wonder whether a sig- 
nificant parity dependence remains in the long-junction 
regime. 

Remarkably enough, the parity dependence becomes 
even stronger. While in a short junction the two branches 
/+((/>) = —/_((/<) differ only in sign, we find that in a 
long junction they differ both in sign and in magnitude. 
In particular, the largest current that can flow without 
dissipation is twice as large for /_ than it is for /_|_. The 
difference is illustrated in Fig. [T] in the zero-temperature 
limit. The basic physics can be explained in simple terms, 
as we will do first, and then we will present a complete 
theory for finite temperature and for arbitrary ratio i/^. 

We set the stage by summarizing the findings of Fu 
and Kane j7j in the short-junction regime. The spec- 
trum of the Bogoliubov-De Gennes Hamiltonian -ffedc 
is a ±e symmetric combination of a discrete spectrum 
for |e| < A and a continuous spectrum for |e| > A. 
Since backscattering along the quantum spin-Hall edge 
is forbidden by time-reversal symmetry [19] . this is a 
ballistic single-channel Josephson junction. In the limit 
L/^ — ^ the discrete spectrum consists of a pair of levels 
at e± = =fA| cos(0/2)|, while the continuous spectrum is 
(/)-independent [50]. Quite generally, an eigenvalue e(0) 
of -ffedG contributes to the supercurrent an amount 
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r£(0), 



(1) 



with g a factor that counts spin and other degeneracies 
|21) . There is no spin degeneracy at the quantum spin- 
Hall edge (since spin is tied to the direction of motion), 
so g = 1 and the level s± contributes a supercurrent [7] 



±tsin(0/2). 



< TT. 



(2) 



To discuss the fermion-parity anomaly we assume, for 
definiteness, that the total number Af of electrons in the 
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FIG. 1: Phase-dependent excitation spectrum of a Josephson junction along a quantum spin-Hall edge (left panels) and 
corresponding zero-temperature supercurrent (right panels). The supercurrent is 47r-periodic, with two branches J+ (blue 
solid), /_ (red solid) distinguished by the ground-state fermion parity and with a parity switch at (fj — ±n. The top row shows 
the short-junction limit of Ref. [71 the bottom row the long-junction limit calculated here. (The jump in /_ at = occurs 
because of the change in slope indicated by the green arrows in the magnified central part of the spectrum.) The 27r-periodic 
supercurrent I2tt without parity constraints is also shown (green dashed). The critical current is the same for and I2tt in 
the short junction, but different by a factor of two in the long junction. 



system is even. (A different choice amounts to a 2vr 
phase shift, or equivalently, an interchange of /_|_ and 
/_.) The ground-state fermion parity a is even for (j) = 
and switches to odd when (j> crosses vr. Since Af is fixed, 
this topological phase transition must be accompanied by 
a switch between even and odd number of quasiparticle 
excitations. At zero temperature only the two levels £± 
closest to the Fermi level (e = 0) play a role, and the par- 
ity switch of a means that a quasiparticle is transferred 
from £+ < to £_ > 0. It cannot relax back from £_ to 
e+ at fixed parity of A/". 

The resulting current-phase relationship can be repre- 
sented by a switch between 2vr-periodic branches I±{<j)) 
(reduced zone scheme), or equivalently as a 47r-periodic 
function l4Tr{4>) (extended zone scheme). Both represen- 
tations are shown in Fig.[l] upper panels. We also include 
the 27r-periodic current /27r that results if the system can 
relax to its lowest energy state without constraints on 
the parity of Af. 

So much for the short-junction limit. An elementary 
discussion of the long-junction regime (to be made rig- 
orous in just a moment) goes as follows. For L ^ ^ we 
may assume a local linear relation between the 

current density I and the phase gradient (f>/L <^ 1/^, of 
the form / = constant x evcj)/ L. The linear increase of /_ 
is interrupted at </> = by a discontinuity A/„ — 2ev/ L. 
Half of it results from the jump in the slope of the low- 
est occupied positive energy level e — {n — \(l)\)hv/2L 
(green arrows in Fig. [l^). The jump in the slope of the 
highest occupied negative energy level contributes the 
other half. In the extended zone scheme, the resulting 



supercurrent I^^-^ is a 47r-periodic sawtooth with a slope 
A/_/47r = eET/2TTh. 

The corresponding parity-dependent supercurrents in 
the reduced zone scheme are 

^+^2^"^' ^- - ^('^-2^sign0), 101 (3) 

The 47r-periodic supercurrent /47r switches from /+ to 
I- sX 4> = vr, while I2tt remains in the branch by 
compensating the switch in ground-state fermion parity 
cr by a switch in the parity of the electron number A/". 
These are the curves plotted in Fig. [l] (1 ower panels). 
The maximal supercurrent is reached near (j) — 2tt for 
(with parity constraint) and near <j) — n for l2-n (with- 
out parity constraint). There is a factor of two difference 
in magnitude of these critical currents in a long junction, 

hrr.c^eET/h, l2^,c^eET/2h. (4) 

In contrast, for a short junction both are the same (equal 
to eA/2H). 

To determine the crossover from the short-junction 
limit ([2| to the long-junction limit ([s]), including the tem- 
perature dependence, we adapt the scattering theory of 
the Josephson effect [25] to include the fermion parity 
constraints. Input is the scattering matrix sq of electrons 
in the normal region and the Andreev reflection matrix 
TA at the normal-superconductor interfaces. These take 
a particularly simple 2x2 form at the quantum spin- 
Hall edge, but our general formulas are applicable also 
to multi-channel topological superconductors. 
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The parity-dependent partition function is [T^HHl [25] 



Z± 



\e>0 / \e>0 e>0 

- iZo('l±ntanh(^e/2)Y 



e>0 



(5) 



with /3 = l/fceT and Zq = ne>o 2 cosh(/3e/2) the par- 
tition function without parity constraints. From the ex- 
pression for Z± one can see that the ± selects terms that 
contain an even (+) or an odd (— ) number of quasiparti- 
cle excitation factors e~^^, as is dictated by the ground- 
state fermion parity. The partition function Z gives the 
free energy F and hence the supercurrcnt / j27| . 



(6) 



/2.^/o = y^, f^o--r'lnZo. (7) 

The density of states p{e) contains both the discrete 
spectrum for |e| < A (a sum of delta functions at the 
Andreev levels) and the continuous spectrum for \e\ > A, 
including also a contribution pg from the superconduct- 
ing electrodes. Scattering theory gives the expression [5S] 



p(e) = Im — i^fe + iO 
de 



-Ps{e), (8) 

= _7r"MnDetX(e), X = (1 - M)Af"i/^ (9) 
Mis) = rl{-e)sl{-s)rAis)so{e). (10) 



The factor M~^/^ in the definition of X, as well as the 
term ps, give a 0- independent additive contribution to 
Fq without any efi'ect on Iq, but we need to retain these 
terms here because they do enter into the parity con- 
straint for 

In the absence of parity constraints, Ref. [551 gives the 
free energy 



^o = -r'E^olnDetX(^Wp), 



(11) 



as a sum over fermionic Matsubara frequencies ujp — 
{2p + l)7r//3. A similar calculation [29 gives the parity- 
dependence in the form 



F^^Fo- In 1 



1 



ae 



^VDetX(O) 



X expl ^(-l)PlnDetX(irjp/2) 
cr = sign [Pf (rASo - Sor1){Detisoy^^^]^^„, 



(12) 
(13) 



with bosonic Matsubara frequencies ftp — 2pn/ ji. The 
ground-state fermion parity a is given in terms of the 
Pfaffian of the anti-symmctrized scattering matrix, eval- 
uated at the Fermi energy. The sign ambiguity in the 
square root is resolved by fixing cr = 1 at = 0. 



Eq. ( 12 ) contains a contribution from the supercon- 



ducting electrodes. 



Js ~ I d£/9s(e) lntanh(/3e/2), 

'A 



(14) 



which only plays a role at temperatures T > A/fce. The 
factor e"'^ can therefore be replaced by unity in the long- 
junction regime, when k^T < i?x ^ A. 

We now specify these general formulas for the quantum 
spin-Hall edge, with Hamiltonian (30) 



HBdG 



vpaz + U [x) 



A" 
vpaz 



{x)ay 
-U{x) 



(15) 



The edge runs along the x-axis, p = —ihdx is the mo- 
mentum operator, and the electrostatic potential is U{x) 
(measured relative to the Fermi level). The pair poten- 
tial A(a;) vanishes in the normal region |a;| < L/2. In 
the two superconducting regions we set A(a:) = Ae±^'^/2, 
with a step at x = ±L/2. This socalled "rigid boundary 
condition" is justified for a single channel coupled to a 
bulk superconducting reservoir IQ' . 

A mode-matching calculation gives the scattering ma- 
trices 



«o = \ „iy n ) ' X(e) = Xo + e/ET, 



e'x 

-ae-^'l'/^ I ' 



(16) 



a{e) 



le 



Det X(e) = 2 cos + a^e"''""^ + a-^e-^^'""' 



A2 ^ A' 

(17) 



We discuss the various terms in these expressions. The 
electron scattering matrix sq is purely off-diagonal, be- 
cause of the absence of backscattering along the quan- 
tum spin-Hall edge. The transmission phase x depends 
linearly on energy because of the linear dispersion. Elec- 
trostatic potential fluctuations contribute only to the 

energy- independent offset xo = ~{^''^)~^ jio U dx, which 
drops out in Eq. (|9|. The Andreev reflection matrix ta 
(from electron to hole) is unitary below the gap. Above 
the gap there is also propagation into the superconduc- 
tor, so TA is sub-unitary. The same expression ( 16 ) for 
ta applies at all energies, evaluated at e + iO+ to avoid 
the branch cut of the square root. 

Putting all pieces together [^S] we obtain the parity- 
dependent supercurrent for arbitary ratio A/Et. In the 
short-junction limit A/Et — >■ we recover the known 
result ([2]), when the energy dependence of the scattering 
matrix and the phase sensitivity of the continuous spec- 
trum can both be ignored. In the opposite long-junction 
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FIG. 2: Phase dependence of the parity-constrained super- 
current (solid curves, in units of eEx/h oc l/L), calcu- 
lated by a numerical evaluation of the Matsubara sums. The 
left panel shows the crossover from the short-junction to the 
long-junction regime in the zero-temperature limit (full in- 
terval —2-71 < (j) < 2tt). The right panel shows the tempera- 
ture dependence in the long-junction limit (reduced interval 
< (f} < 27r). The left panel also shows the supercurrent 
I2tv without parity constraints (dashed curves). The insets in 
the right panel show current-biased superconducting circuits 
that measure the I-V and /-$ relationships of a Josephson 
junction. 

limit A/Et — > oo we find 

= /o - + cos(</)/2)e^--/2^^-] , (18) 

OO 

S = ln(l + 2e-^^/^-^ cos cj) + e-^^"/'^'^) , (19) 

p=i 

2e °° 

^ /q = — sin [cos (/)-!- cosh(2a;p/£'T)]~ . (20) 

^ p=0 

The plot of the results in Fig. [2] shows that the crossover 
from a sine to a sawtooth shape occurs early: already for 
A = Et (so for L = ^) the maximum of the current-phase 
relationship is close to (j) — 2tt. The sawtooth shape is 
preserved with increasing temperature for ksT < \E't. 

These are encouraging results for the experimental ac- 
cessibility of the long-junction regime. The quantum 
spin-Hall effect has been observed in HgTe/CdTe quan- 
tum wells [31], and more recently in InAs/GaSb quantum 
wells [35] — where also Andreev refiection from super- 
conducting Nb electrodes was demonstrated [33]. For 
a typical Fermi velocity of w ~ 10^ m/s in a semicon- 
ductor and superconducting gap A ~ 1 meV in bulk 
Nb, the coherence length is ^ = 70 nm, so the Joseph- 
son junction length L = 0.5 /im from Ref. [33 is deep in 
the long-junction regime. Since the long-junction regime 
is already entered for L w ^, this would apply even 
if the effective superconducting gap is well below the 
bulk value of Nb. The corresponding Thouless energy 
is i^T/fce = 1.5 K, so at T = 100 mK one should be close 
to the low-temperature limit. 

In the ongoing search for the 47r-periodic Josephson ef- 
fect the first results have been reported 34J for the AC 



effect (fractional Shapiro steps [51 [T3HT5] ). A DC mea- 
surement of the current-flux (/-$, 4> = 2e^/K) relation- 
ship, on time scales large compared to the time Tqp ~ /is 
for unpaired quasiparticles to tunnel into the system 35] , 
will measure the 27r-periodic I2tx rather than /4^. Such a 
phase-sensitive measurement (Fig. [2] upper inset) would 
produce the critical current /27r,c without any signature 
of the parity anomaly. In contrast, a p\ids,e- insensitive 
measurement of the critical current through the current- 
voltage {I-V) characteristic (lower inset) will produce 
Iatt,c even on time scales ^ Tqp, because the phase of 
a resistively shunted (overdamped) circuit can adjust to 
a change in M on time scales much smaller than Tqp. A 
change in the parity of A/" will be compensated by a 27r 
phase shift, without a change in critical current [29 . In 
a short junction /2Tr,c and /47r.c are the same, so this does 
not help, but in a long junction they differ by up to a 
factor of two. 

In conclusion, we have presented a theory for the 47r- 
periodic Josephson effect on scales large compared to the 
superconducting coherence length. A multitude of sub- 
gap states, as well as a continuum of states above the 
gap, contribute to the supercurrent for L ^ ^, but still 
the parity anomaly responsible for the 47r-periodicity per- 
sists. In fact, we have found that in a long junction the 
anomaly manifests itself also in a phase-insensitive way, 
through a doubling of the critical current. This opens up 
new possibilities for the detection of this topological ef- 
fect at the quantum spin-Hall edge f3P-'33] , and possibly 
also in semiconductor nanowires [34, 36-38 . 
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knowledged. This research was supported by the Dutch 
Science Foundation NWO/FOM and by an ERG Ad- 
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Appendix A: Details of the calculation of the free 
energy 

1. Transformation from the real to the imaginary 
energy axis 

According to Eq. ^ the free energy Fa- — — inZ^, 
with cr = ±1 the ground-state fermion parity, is given by 



£>0 



1 + a exp ^lntanh(/3e/2) 
Fo = -/?-^ ln[2 cosh(^e/2)] 



(Al) 
(A2) 



e>0 



Here Fq is the free energy in the absence of parity con- 
straints. The infinite product over e is defined in terms 
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of the density of states p by 



exp 



E>0 



dep{e)\TLf{e) 



(A3) 



with ojp = (2p + l)Tr//3 the fermionic Matsubara fre- 
quency. For Fa- the Matsubara sum contains bosonic 
Matsubara frequencies flp — 2pn/l3, 



Integrals of the type (A3) can be done efficiently by 
contour integration, made possible by the fact that the 
scattering matrices are analytic in the upper half of the 
complex energy plane. This was worked out in Ref. |28] 
for Fq, and here we adapt that method to include the 
parity constraint in F^. 

The function ly in the density of states (|8| satisfies 
p{e) = v*{—e), expressing the electron- hole symmetry. If 
f{e) is an even function of e, we may therefore convert the 
sum X]e>o /(^) over positive energies (including both the 
discrete and the continuous spectrum) into an integral 
along the entire real energy axis of fdu/de, or vdf /de 
after a partial integration. Closing the contour in the 
upper half of the complex energy plane picks up the poles 
of df /de, which for the free energy are the Matsubara 
frequencies on the imaginary axis. 

In Ref. 28 this transformation was carried out for Fq, 
leading to 

/•CO 

Vln[2cosh(/3£/2)] = / de p(e) ln[2 cosh(/3e/2)] 
dv{£ + iO+) 



6>0 



1 

2i 



de 



— oo 
oo 



de 



■ln[2cosh(/3e/2)] 



dei/(£ + i0+)tanh(/3e/2) 



1 

^lntanh(/3e/2) = tt. i de 



di^{e + iO+) 



de 

X ln|tanh(/3e/2)| 



I 
2i 



dev{e + «0+)- 



1 



sinh /?£ 



-TT Y^{-lYv[inp/2) - \ttv{0). (A5) 
p=i 



The notation j indicates the Cauchy principal value of 
the integral (which picks up half of the pole at £ = 0). 
We have made use of the identity 



/ 



dx -f- In \x\ 
dx 



j-dxfix)^. 



(A6) 



-TT^l'iiujp), 

p=0 



(A4) 



For the free energy this gives 



Fa = -^"i^lnDetX(iWp), 

p=0 



Fa^Fo- In i 



H-cre'^VDetX(0)exp ^(-1)^ In Det X(iOp/2) 
\p=i 



(A7) 
(A8) 



where we have substituted i/ = — 7r~^lnDetX from 
Eq. ([9| and also included the factor e''^ from the cj)- 
independent density of states in the superconducting 
electrodes. 



2. Scattering formula for the ground-state fermion 
parity 



(T-\/Dct X{0). We seek to express this quantity in terms 
of the normal-state scattering matrix so(0). We will 
then have obtained a scattering matrix expression for the 
topological quantum number a — equivalent to Kitaev's 
Hamiltonian expression [B] but requiring only Fermi-level 
information. 



The free energy (AS) incorporates the ground- 
state fermion parity dependence through the quantity 



We start from the definition of X = (1 - M)M-i/2 
and use that M(0) = U*U with unitary U = rA(0)so(0). 
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Since Det M = 1 , we have 

DetX = Det(l - U*U) = Det (1 - U*U)'^ 
_ Det jU - U^) _ [Pf ([/ - 
~ DetU ~ BeiU 
The determinant of U is independent of (f>, 

Det [/ = Det [rA(0)so(0)] = Det [iso(O)]. 



(A9) 



(AlO) 



We may now identify 



a^/'DetX{0) 



sign 



_ PfjU^-U) 
Pf(rASo 



\/Det isQ 



(All) 



The two branches of the square root function introduce a 
sign ambiguity, which is resolved as follows. The square 
root of Det X > is taken on the principal branch, while 
the branch of the square root of Det iso is fixed by set- 
ting fT = 1 at = 0. Once this sign is fixed, the phase 
dependence of the topological quantum number i7((/)) is 
determined by Eq. (All ) entirely in terms of Fermi-level 
properties. 



The scattering formula (All) is applied in the main 
text to a single-channel Josephson junction, where we 
know a priori that the topological phase transition oc- 
curs at = TT (mod 27r). More generally, in the multi- 
channel case there may appear any even number of sign 
changes of a between = and (j) = tt. Eq. (All I is then 



particularly useful, since it permits the determination of 
a at a given (j) without having to track the number of 

sign changes back to (/> = 0. 

We conclude this subsection by placing Eq. ( All[) in 
the general context of topological states of matter |39]. 
The ground-state fermion parity a is the Z2 topological 
quantum number of a system of dimensionality d = 
in symmetry class D (when sq has no symmetry restric- 
tions) or BDI (when sq = sj). The dimensionality zero 
refers to the fact that this is a closed system. We may 
consider opening up the system, promoting it to d = 1, 
by replacing one of the two superconducting contacts 
by a normal metal. The topological quantum number 
Q then counts the number of Majorana zero-modes at 
the normal-superconducting interface. It is given by the 
determinant of the reflection matrix (for class D, with 
Q G Z2) [40] or by its trace (for class BDI, with Q € Z) 

m- 



3. Regularization 



The transformation of the integral ( A5 1 over real ener- 



gies into a sum over imaginary frequencies requires that 
v{iuj) goes to zero faster than l/w for ui — >■ 00. To ensure 
this, we decompose v = j/qo + ^v^ with v^o the large-w 
limit of u[iuj). It is convenient to specify i^oo(O) = 0. The 
integral over Vao is done along the real energy axis, where 



it converges, and then the remaining integral over 5v be- 
comes a converging sum over Matsubara frequencies. 

More specifically, for the quantum spin-Hall edge we 
take 



,{e) ^ -^-Mn[(l-4eVA')e- 



-2ie/ET^ 



(A12) 



in view of Eq. (17 1. The integral over Voo can be evalu- 



ated in closed form, 

J DO 



2i 



d£i'oo(£ + iO+)— — 
smh pe 



1 



TT 



dx 



sinha; 



dx ■ 



1 



/3A/2 



sinhx 



lntanh(/3A/4). 



(A13) 



The remaining integral over 5u = v — Vao then becomes 
a convergent sum over Matsubara frequencies. 



2i 



sinh/3. =^-1^0). 



5 = -7r^(-l)P[iy(zr!p/2) - z.oo(*^^p/2)]. 

p=i 



(A14) 
(A15) 



This gives the regularized version of Eq. ( A8 ) 



= Fq - 13-^ In i (1 + (76'^+"'°°+'^^ ^/BetX{0) ) . (A16) 



4. Evaluation of the supercurrent along the 
quantum spin-Hall edge 



We apply Eq. (AI6) to the quantum spin-Hall edge. 
According to Eq. (17), we have 



^{iLj) - 

C±(w) 



1 



ln[2cos0-f C+l'^) + C-(w)], (A17) 



(v/l+a;VA2±w/A)^ 



(A18) 



The Matsubara sum (A15l, with given by Eq. ( A12[ ), 
takes the form 



5 = ^(-l)nn 



2 cos C+ (^^p/2) + C- (^^p/2) 



(1 + 1]2/A2)ef2p/^^T 



The function Jrx, is given by Eq. (A13) and 



(7^T)ciX{Q) = 2 cos(0/2). 



(A19) 



(A20) 



in view of Eqs. (16) and ( |All| ). 

The superconducting electrodes couple to the quantum 
spin-Hall edge via a single transverse mode, over a total 
length Ls ■ The corresponding density of states is 



Ps(£) 



t:Es - A2 ' 



\e\ > A. 



(A21) 
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We have defined Eg = hvp/Ls- The supercond ucting 
electrodes affect the parity-dependent free energy (A16) 
through the factor e"^^, with 



de ps(e) lntanh(/3e/2) 



2/3 



TrE< 



de 



Ve 



A2 



S J A 



sinh Pe 



(A22) 



Collecting results, we arrive at the parity-dependent 
supercurrent 



<7 — sign [cos(0/2)], 



(A23) 
(A24) 



with Iq the 27r-pcriodic supercurrent in the absence of 
parity constraints, 

/o = T^sin</)^[2cos0 + C+K)+C-K)]"'- (A25) 

In the long-junction regime A ^ E'y, k^T these results 
reduce to the equations ( 18 )-( 20 1 given in the main text. 



5. Short-junction limit 

As a check on the consistency of the whole formalism, 
we take the short-junction limit Ej^ — > oo of the parity- 
dependent supercurrent (|A23 ) and see if we recover the 
results of Fu and Kane [7j. We choose the phase interval 
loil < TT and abbreviate 



u= \cos{(j)/2)p/\. 



(A26) 



The Matsubara sums (A19) and (A25) can be evalu- 



ated in closed form in the short-junction limit. 



lim 5 = V(-l)Pln 



Et— ^OO 



lim Jo 

Et— >-oo 



2 cos ( 



tanh '( 



1 + r!2/A2 



\ 2 cos((/>/2) tanh(/3A/4) ) ' ''^ 

4e 



hp 

eA 
2S 



sm (f) [2 cos (/) - 

p=0 

sin(0/2) tanhu. 



(A28) 



In the same limit J^o ~ lntanh(/?A/4); upon substitution 
into Eq. ( A23 1 we arrive at 



/,=/o-||ln(i±ie^stanh, 



2e d 
hPdcj) 



ln(coshu ± sinh 



(A29) 



In the zero-temperature limit Js — ^ and we recover 
the result of Ref . 

lim /± = T ^ - ± ^ sin((/./2) . (A30) 

The parity dependence at finite temperature can be 
quantified by the difference 61 — |(/+ — /-), for which 
we find the compact expression 

in agreement with Ref. 11^ 



6. Zero-temperature limit in the long-junction 
regime 

Another check on the formalism is provided by the 
combined zero-temperature and long-junction limits. We 
again choose the interval < tt. The Matsubara sums 
(A19) and (A25) are given in the long-junction limit by 



lim 5 = ^(-l)Pln(l + 2e^"''/-^Tcos< 



-2^2p/_ET^ 



A-s-oo 



p=l 



lim Iq = ^ sin (f> [cos (j) + cosh(2ajp/i?T)] ^ ■ 

A-i-oo hp ^-^ ■' 



(A32) 



(A33) 



In the zero-temperature limit the sums may be con- 
verted into integrals, with the results 



5 ^ -ln|2cos(0/2)|, /q 



2TTh 



(A34) 



The two terms Joo and Js both vanish at T = 0. Sub- 
stitution into Eq. (A23) gives /+ = /q, in agreement 



with Eq. ([3]), while /_ remains undetermined. The zero- 
temperature limit of /_ depends on on higher order 
terms in the low-temperature expansion of S that we 
have not been able to calculate analytically. A numer- 
ical calculation (using the Pade approximant built into 
the NSumI AlternatingSigns routine of Mathematica) 
gives 

(A35) 

resulting in a current /_ in agreement with Eq. (|3|. 



Appendix B: Circuits to measure the critical current 
with and without parity constraints 

As explained in the main text, the two circuits shown 
m Fig. [2] (inset) both measure the critical current of 
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a) 




FIG. 3: Left panel a: Circuit of a current-biased, resistively 
shunted Josephson junction, to measure the current-voltage 
characteristic. Right panel b: Circuit of an rf SQUID to mea- 
sure the current-phase relationship. 



Et. We also wish to ensure that the Josephson junction 
remains in its ground state during the phase relaxation, 
which requires Ej^Tj/h » 1. Together these three condi- 
tions are met if 

{R/RqfTj < RC « Tj, (B3) 

with i? < i?Q. 

At a fixed parity V, the bias current I = I-p + V/R 
drives the phase 4>{t) according to 



the Josephson junction, but their sensitivity to fermion- 
parity constraints is fundamentally different. A measure- 
ment of the current-phase relationship (upper circuit) is 
insensitive to parity constraints when quasiparticles can 
enter or leave the system on the time scale of the mea- 
surement. This gives the critical current I2tt.c- A mea- 
surement of the current-voltage characteristic (lower cir- 
cuit) remains governed by fermion- parity constraints as 
long as the quasiparticle tunneling time Tqp is large com- 
pared to the phase relaxation time tj of the resistively 
shunted Josephson junction. This then gives /47r,c- Here 
we analyze these two circuits in some more detail. 

In the zero-temperature, long-junction limit, we 
have the 47r-periodic sawtooth current-phase relationship 
shown in Fig. [l] (lower panel) . This plot is for an even 
number of electrons in the system, V = (—1)'^ = 1, while 
for odd parity V = —1 the sawtooth is displaced horizon- 
tally by 2tt. Both cases are contained in the formula 



■mod4Tr(0 + Vtt + tt) — I^TT,, 



(Bl) 



4?™, 



The modulo function is defined by mod47r('/') = 
with n e Z such that mod47r((^) € [0,47r). 

We start by considering the current-biased circuit of 
Fig. [3^. A voltage V = {h/2e)d(f)/dt drops over a resistor 
R in parallel with the Josephson junction, of capacitance 
C. The two characteristic time scales of the circuit are 
the RC time and the phase relaxation time 



h 1 
2e RliTT I 



Rq±_ 
2R E'-£ 



(B2) 



where Rq — h/e^ is the resistance quantum. The char- 
acteristic energy scales of the Josephson junction are 
the charging energy e^/C and the Josephson energy 
fih-K.c/e = Et. 

The capacitance should be suffiently small that the 
phase dynamics is overdamped, RC <^ tj, and suffi- 
ciently large that the phase dynamics is classical, e^/C ^ 



I Tj d(f> 



2tt 



-mod. 



Vtt 



(B4) 



A typical value E-r/kB ^ 1 K gives H/Et ^ 10"" s. The 
typical time scales for quasiparticle poisoning are in the 
^s to ms range [3S], so even if R <C_Rq we can safely 
assume that tj ^ Tqp and use Eq. (B4) to calculate the 



relaxation of the phase in between two quasiparticle tun- 
neling events. 

The phase relaxation due to a quasiparticle tunneling 
event at i = (by which V i-> —V) amounts to a 2tt 
phase slip on a time scale tj, 

cp{t) = (/.(O)e-*/"'' + [0(0) + 27r](l - e"*/"''). (B5) 

Before and after the phase slip the junction is in a zero- 
voltage state, for bias currents / < l4n,c- During the 
phase slip there is a voltage pulse of integrated area 
/ V{t)dt = h/2e. The corresponding time-averaged volt- 
age V = h/2eTqp is smaller by a factor Tj/rqp <^ 1 than 
the voltage that develops for / > lii^^c- 

This shows that the current-biased circuit of Fig. [3^ 
provides a DC measurement of the parity-constrained 
critical current I/t-K^c- In contrast, the circuit of Fig. 3j3 is 
not sensitive to parity constraints. This rf SQUID is phase- 
biased for sufficiently small inductance L <^ h/el4Tr,c- 
Quasiparticle tunneling events have only a small effect 
on the phase, which remains fixed by the enclosed flux 
$ = (fi/2e)(/) « LI. 

At low temperatures the parity of the number of elec- 
trons J\f in the system will equilibrate at the ground- 
state fermion parity a, which implies that Af is even 
{V ^ 1) for mod4^((/) -f tt) < 27r and odd {V = -1) 
for mod47r(0 + n) > 2tt. In either case the supercur- 



rent I-p given by Eq. (Bl| cannot become larger than 
lA-n,c/2 = /27r,c — which is the critical current without 
parity constraints. 



